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ABSTRACT 

Simulations of planetary microlensing at high magnification that were carried out on a 
cluster computer are presented. It was found that the perturbations due to two-thirds 
of all planets occur in the time interval [— O.StpwHM, 0.5<fwhm] with respect to the 
peak of the microlensing light curve, where ipwHM is typically ~ 14 hours. This implies 
that only this restricted portion of the light curve need be intensively monitored for 
planets, a very significant practical advantage. Nearly all planetary detections in high 
magnification events will not involve caustic crossings. We discuss the issues involved in 
determining the planetary parameters in high magnification events. Earth mass planets 
may be detected with 1-m class telescopes if their projected orbital radii lie within 
about 1.5 — 2.5 AU. Giant planets are detectable over a much larger region. For multi- 
planet systems the perturbations due to individual planets can be separated under 
certain conditions. The size of the source star needs to be determined independently, 
but the presence of spots on the source star is likely to be negligible, as is the effect 
of planetary motion during an event. 

Key words: Gravitational lensing: microlensing - Methods: numerical - Stars: plan- 
etary systems 



1 INTRODUCTION 

New results on extra-solar planets are now being obtained 
by various means, most notably the radial velocity, tran- 
sit and gravitational microlensing techniques (Ferryman, 
2000). Gravitational microlensing events of high magnifica- 
tion have been shown to be promising in this regard (Griest 
& Safizadeh, 1998; Rhie et al, 2000; Gaudi et al, 2002; Bond 
et al, 2002a). Using new observational and analytical tech- 
niques, it has recently been shown that high magnification 
events may be readily detected, including events with very 
high magnifications (Bond et al, 2002b). Here we discuss 
the detectability of planets in this new class of events using 
1-m class telescopes, including: 

• Critical observational period for planet detection 

• Zones of detectability of terrestrial and giant planets 

• Detectability of habitable planets and solar-system ana- 
logues. 

• Precision of planetary system characterisation 

• Multi-planet systems 

• Effect of the source star radius and stellar spots 

• Effect of orbital motion during a microlensing event 

• Comparison with low magnification events 

• Absolute characterisation of planetary parameters 



2 DEFINITIONS USED IN THIS STUDY 

We consider the lensing of a source star by a lens star where 
the lens star is assumed to have a planet or planets orbiting 
it not dissimilar to those in our solar system. To first order, 
the lensing by such a system can be approximated by the 
lens star alone. In this case the lensing can be parametrised 
in terms of the Einstein radius Re and crossing time 
where 



Re = 4.42 



and 



tE = 38.25 



Ml 



O.3A/0 V 8kpc 



Dos 



^J(l-d)d AU 



Ml 



O.SMq V 8kpc 



Dos 



(1 ~ d)d days. 



(1) 



(2) 



-DoL and Dls are the observer-lens and lens-source dis- 
tances respectively and d — . We see that Re — 1.9 
AU and tE — 16.6 days for bulge events when the lens mass 
Ml — O.3M0, -DoL — 6kpc, Dos — 8 kpc. The transverse 
velocity, vt, of the lens with respect to the observer- source 
line of sight is taken as 200 kms"'^. 

High magnification occurs when the lens and source 
stars are well aligned, i.e. when the impact parameter 
■"min ^ -Re- In this case 



A, 



Re 



> 1 



(3) 
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A useful quantity in high magnification events is the time of 
full-width-half-maximum with respect to the peak given by: 

3.5fE 



tpWHM 



(4) 




e = -30 



3 GENERAL STRUCTURE AND CRITICAL 
TIME OF PLANETARY DEVIATIONS IN 
HIGH MAGNIFICATION EVENTS 

Liebes (1964) discussed several features of gravitational mi- 
crolensing events with magnifications up to ~ 1000, and 
remarked that planets orbiting the lens stars in these events 
would perturb the light curves. Griest & Safizadeh (1998) 
subsequently computed the perturbations for a variety of 
high magnification configurations, and showed that the 
probability for detecting planets in these events is high. In 
particular, they showed that planets with masses as low 
as IOMq could be detected with significant probability in 
events with magnifications ~ 50 by monitoring the peaks of 
the events with a photometric precision ~ 1%. 

It has subsequently been shown that events with mag- 
nification up to ~ 200 are quite readily detectable, and also 
that intensive photometric monitoring of the peaks of these 
events with a precision significantly better than 1% is feasi- 
ble (Bond et al, 2002a; Bond et al, 2002b). Consequently 
the prospects for planet detection have improved. 

In this section we examine typical light curves for events 
of higher magnification, to determine the basic requirements 
for planet detection. The light curves presented here, and 
indeed throughout the paper, were generated using a rapid 
version of the inverse ray shooting method (Schneider & 
Weifi, 1986; Kayser, Refsdal & Stabell, 1986; Wambsganss, 
1997). The technique was implemented on a ~ 200 node 
cluster computer at the University of Auckland (Dobcsanyi, 
1999). A further improvement in the code runtime is de- 
scribed in Appendix A. 

In Figure llj we show a selection of light curves for a typ- 
ical lens system consisting of a star with a single planet in a 
typical event of high magnification. The fractional difference 

(a —A- * 

light curves are plotted: 5 = 



where 



..ii^i^ ^~Lp is the light 

curve of the lens star including its planet and As is the light 
curve of the lens star alone. The planet-star mass ratio of 
the lens is e = 4 x 10"''. This corresponds to a planet of four 
Earth masses orbiting a O.SMq mass star, the most probable 
lens mass. The maximum magnification is j4max = 80. The 
radius of the source star, Rs, has been set equal to Rq, as 
is done throughout this paper unless otherwise stated. The 
projected orbital radius of the planet at the time of the lens- 
ing event is 0.8i?E or ~ 1.5AU. The lens geometry is shown 
in Figure ^ 

It is clear from Figure |l| that the position angle of a 
planet at the time of lensing may be determined from the 
light curve. This follows because the peak deviation occurs 
at the point on the source star track when it is crossed by 
the planet-star axis. Thus the position angle 9 is given by 
the equation 



tan( 



VTtp 



(5) 
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where tp denotes the time of maximum planetary devia- 



tion. 



Equation ^ implies that, provided 6 is not within 



Figure 1. Sample perturbed light curves for a single planet with 
mass ratio e = 4 X 10~^. Each graph shows the deviation from 
the single lens light curve. The maximum light curve amplifica- 
tion is Amax = 80. The angle between the source star track and 
the planet-star axis is denoted by 9 and defined in Figure bl For 
negative values of d the source star track crosses the planet-star 
axis between the planet and star. The projected planet orbit ra- 
dius is 0.8-Re- For 9 = 0° and 180°, there is a planetary caus- 
tic crossing, resulting in a large deviation from the single lens 
light curve. The time axis is in normalised time units, % = t^. 
In these units tpwHM = 0.044 when Amax = 80. The interval 
[— 0.5tpwHM, O.StpwHM] is indicated on the lower-most set of x- 
axes as a thick line. 



±30° of either 0° or 180°, then tp lies in the interval 
[— O.StpwHM, 0.5tFWHM]. It follows that for two-thirds of all 
planetary positions, the perturbations they produce to the 
microlensing light curve occur in the very restricted time in- 
terval [— O.StpwHM, O.StpwHM]. This is the critical time for 
planet searches. 

Figure ^ shows typical examples of deviations produced 
by planets of different masses and at different projected 
radii. In general, heavier planets produce larger deviations, 
as do planets nearer the Einstein ring. The shapes of the 
deviations are also dependent on these properties. It is thus 
possible to determine both the mass and the projected ra- 
dius of a planet from the deviation to the light curve. In 
Section ^ we discuss the precision with which this may be 
achieved. 
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Figure 2. Co-ordinate system used by the inverse ray shooting 
code. The scales are in units of the Einstein ring radius, i?E- 
In what follows the projected orbit radius at the peak of the 
microlensing, y^x'^ + (j/p — Wmin)^i is denoted by a. 
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Figure 3. Light curve perturbations for various planet-star mass 
ratios and separations. The maximum magnification A^ax = 100, 
and e = 60°. 



4 PLANET DETECTABILITY 

Planetary zones of detectability in high magnification events 
were determined by comparing simulated perturbations to 
microlensing light curves produced by planets of various 
mass and at various projected radii with minimum de- 
tectable perturbations. 



4.1 Simulation details 

Each simulated planetary light curve was comprised of 301 
evenly spaced points in the critical time interval Ifwum- In 
order to simulate the observation process, noise was added 



to the simulated light curves. As a representative example, 
we assumed observations of 200s duration were made with 
a Im telescope in a passband of 200nm. In this case, the 
number of detected photons per exposure is given by the 
expression (Eccles 1983): 



photons 



(6) 



where / denotes the telescope throughput. The number at 
any magnification, given a maximum magnification, A^^ax , 
is: 

A 



N{A) = Nn 



Art 



photons. 



(7) 



If we assume Poisson statistics, the photon noise, aN, is 
[N{A)]~^. Using Equations ^ and this becomes: 



aN{A) 
where m. 



nf X 10' 



(-0.4m„,ax-l-ll-5) 



An 



(8) 



'max is the magnitude at peak magnification. Equa- 
tion ^ defines our noise function for any light curve with 
peak amplification Ama.^ and peak magnitude mmax. The 
minimum error that the difference imaging technique can 
achieve is typically: 



(A) = V2a^{A) 



(9) 



because it involves the difference of two similar quantities. 
In what follows we assumed / = 0.3 and Gaussian noise with 
crfj = 2aN. This implies ct^ = 0.36% at the peak of a light 
curve with mmax ~ 15, consistent with the results reported 
by Bond (2002a). 



4.2 Results 

For each light curve was computed using the single lens 
light curve. This was compared to Xoi value obtained 
from data without any planetary signal. A deviation in the 
light curve was considered detectable if the difference be- 
tween these two values, Ax^ = X'^ ~ Xoi greater than 
60. Figure ^ shows the positions for an Earth-mass planet 
which give Ax^ > 60. The threshold value of 60 was chosen 
because it corresponds to ~ 20 consecutive deviations of 1- 
2a (Bond et al, 2002a). For events with ^max ~ 100, the 
total number of measurements that can be made in the time 
interval [— O.SipwHM, O.SfpwHM] is ~ 300 assuming parame- 
ters as above and a readout time between exposures of the 
CCD camera ~ 100 sees. If measurements are made over 
the interval [— tpwHM, tpwHivi], then ~ 600 measurements 
may be anticipated. The probability for a random fluctua- 
tion to yield Ax^ > 60 is < 1% for 300 measurements, and 
~ 4% for 600 measurements. However, any random fluc- 
tuation will be easily identiflable in any densely sampled 
light curve, because it will not produce a coherent pertur- 
bation. The detection limits for the microlensing situations 
that were tested are shown in Table |l]. 

The trends in the results are clear. The detection limits 
extend around the Einstein radius in the lens plane with 
heavier planets detectable further away from this region. 



4.3 Habitable planets 

Numerical simulations by Wetherill 
5 — 15% of all stars with mass 0.5 - 



(1996) predict that 
1.5M0 have a hab- 
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Table 1. Planet-star orbit detection limits. The planet mass ratios used were: 1 X 10-5, 1.7 X 10-* and 3.2 x 10-^. These values 
correspond to Earth, Neptune and Jupiter mass planets, given a O.SMq lens star. Each simulation light curve was comprised of 301 
points over the interval: [— ^tpwHMi f^FWHlvl]- Two values of maximum light curve magnitude are used, and three values of maximum 
light curve amplification. For a perturbation to be considered detectable, the difference in between the perturbed light curve and a 
single lens light curve must exceed 60. The detection limits are in units of Ke- 

Ml = O.3M0 iJs = Rq 

Amax 50 100 200 

mniax 15 17 15 17 15 17 

Earth 0.88 - 1.13 0.95 - 1.05 0.83 - 1.2 0.92 - 1.07 0.8 - 1.3 0.9 - 1.1 

Neptune 0.5 - 2.0 0.7 - 1.6 0.4 - 0.6 - 1.9 0.1 - 0.4 - 

Jupiter 0.0- 0.1- 0.0- 0.0- 0.0- 0.0- 



A/l = O.3M0 = 2R0 



50 100 200 



mmax 15 17 15 17 15 17 

Earth 0.9 - 1.1 0.99 - 1.01 0.87 - 1.15 0.99 - 1.01 0.84 - 1.2 0.99 - 1.01 

Neptune 0.5 - 2.0 0.6 - 1.5 0.4 - 0.5 - 2.0 0.1 - 0.4 - 

Jupiter 0.0- 0.1- 0.0- 0.0- 0.0- 0.0- 



Ml = O.3M0 = 4i?0 



50 100 200 



mmax 15 17 15 17 15 17 

Earth 0.94 - 1.06 - 0.94 - 1.06 - 0.99 - 1.01 

Neptune 0.5 - 0.65 - 1.54 0.3 - 0.5 - 2.0 0.2 - 0.3 - 

Jupiter 0.0- 0.1- 0.0- 0.0- 0.0- 0.0- 



Ml = A/0 = Rq 



Amax 50 100 200 



^max 


15 


17 


15 


17 


15 


17 


Earth 


0.93 - 1.08 


0.94 - 1.06 


0.91 - 1.1 


0.92 - 1.09 


0.88 - 1.14 


0.89 - 1.12 


Neptune 


0.6 - 1.6 


0.75 - 1.3 


0.5 - 2.0 


0.5 - 2.0 


0.4 - 


0.4 - 


Jupiter 


0.0 - 


0.0 - 


0.0 - 


0.0 - 


0.0 - 


0.0 - 



Table 2. Habitable Zone boundaries during the first 10 years 
(Wetherill, 1996; Kasting, 1993). 



Stellar Mass (Mq) 
0.5 
1.0 
1.5 



HZ boundaries (AU) 
0.20 - 0.40 
0.82 - 1.40 
1.70 - 2.80 



itable planet. We carried out simulations to determine if 
this interesting prediction could be tested in microlensing 
events of high magnification. Following Kasting (1993) and 
Wetherill (1996), and allowing for internal heat sources such 
as long-lived radionuclides (Ferryman, 2000), we assumed - 
see Table |^ - the following limits for the habitable zone 
around a star of approximately solar mass: 0.8 — 2 AU. Fol- 
lowing Deeg (2000), the mass range of habitable planets was 
taken to be 0.8 — SM®, although a broader range could have 
been assumed (Ferryman, 2000). 

Simulated hght curves for planets over a wide range 
of masses and orbital radii were generated. Each simulated 
light curve was compared with its associated singe lens hght 
curve, and the maximum deviation, (5max, was calculated. 
These maximum deviations were plotted on the mass-orbit 
plane, see Figure B. Contours for (5max = 0.1%, 0.5% and 



1.0% are shown. In practice, the minimum detectable value 
for 5max for a 1-m class telescope will lie somewhere between 
the 0.1% and 0.5% contours, with the lower limit being more 
closely achieved by a space-borne telescope. We thus con- 
clude that, whilst terrestrial planets are certainly detectable 
in microlensing events of high magnification, habitable plan- 
ets are only marginally detectable, and Earth-Sun analogues 
are probably not detectable using this technique with Im 
class telescopes. 



5 PRECISION OF PLANETARY SYSTEM 
CHARACTERISATION 

To compute the precision at which planetary masses and 
orbital radii may be measured, a single hght curve with a 
planetary deviation was compared to a set of model hght 
curves where the planet position and mass were varied over 
suitable ranges for a typical high magnification event with 
^max ~ 100. As in previous analyses, a degree of merit 
was computed for each comparison curve, assuming mea- 
surements were made at 5 minute intervals, as above. A 
typical plot of the values over the mass-position plane for 
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Figure 4. Detection limit maps for an Earth mass planet orbiting a 0.3Mq lens star. The axes are in units of i?E- The I band magnitude 
at maximum is 15, and the maximum amplification is 100. The map on the left was generated using observation light curves comprised 
of 301 points over the interval [— ^tpwHMi 5*fwhm]- The right map was created using light curves comprised of 601 points over the 
interval [— tpwHMi tpwHlvi]- The detection limits do not change greatly with increased observation interval, but the detection efficiency 
increases slightly. 
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Figure 5. Detection limit contours with the HZ (see text). The 
parameter positions used in the simulations of both an Earth and 
a heavy, cool Earth are denoted by the © and (gi symbols respec- 
tively The parameter values of other solar system planets are in- 
dicated also. The mass ratio and Einstein ring radii are scaled for 
a solar mass lens star. The simulated light curves were computed 
over the time interval [— 2tpwHM; SipwHlvi]) with Amax = 100. 



a planet at (a, e) — (0.8i?E, 5 x 10"^) is shown in Figure | 
Two minima are visible, one at (a, e) = (0.8i?E, 5 x 10" 
which corresponds to the light curve under analysis itself, 
and one at the same mass value, but with a — 1.25i?E. This 
is the well known degeneracy in planet-star distance, where 
a planet at orbital distance a produces a light curve pertur- 
bation which is indistinguishable from that with a planet at 
distance 1/a (Griest & Safizadeh, 1998). There are no other 
minima seen in the mass-position plane, indicating that 
a low mass planet near the Re cannot be modelled by a 
higher mass planet further away. Further simulations using 
other planet mass-position values show that the confidence 
intervals for planet mass and position increase for a planet 
closer to the Einstein ring. Typically, the precision of mass 
measurements for terrestrial planets will be ±50%, and their 
instantaneous projected orbital radii will be determined to 
an accuracy of typically ±5%. The latter measurement will 
be subject to further uncertainty due to the degeneracy men- 
tioned above. 



6 MULTI-PLANETARY SYSTEMS 

Gaudi et al. (1998) commented on the need to consider the 
effects of multi-planet lens systems in high magnification 
events. However, even an initial search for a single planet 
signal in a light curve involves a large amount of computa- 
tion. An initial search for a two planet system which varies 
the minimum number of parameters would require 0(10^^^) 
simulations. Clearly, searching the parameter spaces for a 
multi-planet system is impossible without a priori informa- 
tion. 

Rattenbury (2001) and Bond et al. (2002a) suggested 
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Figure 6. computed on the planet distance-mass plane for a 
planet at (a, e) = (O.SRe, 5x10~^). The angle that the planet-star 
axis makes with the source star track is 30 degrees and Amax = 
100. The contours indicated in the figure are the 1, 2 and 3cr 
confidence intervals, corresponding to Ax^ = 2.3, 6.17 and 11.8 
respectively. 



that if an observed light curve shows evidence of pertur- 
bations due to a multi-planet lens system, it may not be 
necessary to search the planet position and mass parame- 
ter spaces completely in order to identify the lens system 
components. Rather, the more complicated lens system can 
be approximated as a combination of single planet systems. 
A similar result was subsequently reported by Han et al. 
(2001) for low magnification events. Here we investigate un- 
der what circumstances the result is valid in high magnifi- 
cation events. 

A series of two planet system light curves was generated, 
varying the angular separation between the two planets. The 
two planet masses were equal. We then combined the two 
light curves obtained using the corresponding single planet 
systems. The two planet system light curves and the light 
curves from synthesising single planet light curves were then 
compared. It was found that there was little difference be- 
tween these light curves. Figure |^ shows two single planet 
curves, ^i, S2', the actual two-planet light curve and the 
synthetic light curve S = Si + S2. 

A second two planet system was then modelled by a 
summation of the individual single planet light curves, see 
Figure ^. In this case the two planet masses were unequal, 
and a similar result was obtained. A three planet model was 
also tested in the same manner, with similar results. The 
above investigations confirm that multi-planet systems may 
indeed be approximated by a combination of single planet 
systems as a first approximation for more accurate mod- 
elling. 

A further consideration in multi-planet systems is the 
limits that can be set on the resolution of separate planetary 
signals. By taking a series of two planet light curves with 
varying angular separation between the planets, we found 
that the separate perturbations can typically be resolved 
when the planets are separated by > 20 degrees. 



Figure 7. S light curves showing that the synthesis of two sin- 
gle planet light curves is an adequate first approximation to a 
two planet system. On this scale, the "two-planet" and "summa- 
tion" graphs are indistinguishable. In this example, the planets 
are of equal mass (ei,2 = 4 X 10~^), each at 0.8-Re from the lens 
star with a separation of 40° . It is unlikely that two planets will 
have equal, or similar orbit radii, however they can have equal 
projected radii. 
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Figure 8. Synthetic (summation of single planet light curves) 
and multi-planet light curves for a two planet system. The planet 
parameters used are ei = 1 X 10~^, £2 = 3 X 10~^; ai = 0.9Re 
and a2 = 4iJE ■ These parameters correspond to Earth and Jupiter 
mass planets at projected radii 1.7 and 7.6 AU respectively, as- 
suming a O.SA^Q mass lens star. Amax = 100. The planet positions 
are indicated in the figure insert; 9i = 45° and 62 = 135° . 



7 SOURCE STAR EFFECTS 
7.1 Source size effects 

The effect of the size of the source star in planetary mi- 
crolensing was first investigated by Bennett & Rhie (1996). 
These authors considered events of low magnification, and 
found a significant effect. Subsequently Griest & Safizadeh 
(1998), Rhie et al. (2000) and Gaudi et al. (2002) noted 
that a similar effect is present in events of high magnifi- 
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Figure 9. Effect of source star radius on a planetary perturba- 
tion. The units given are in solar radii, assuming a O.SMq lens 
star. The planet parameters are (a, e) = (0.8iJBi4 X 10~^) with 
Amax = 100 and 9 = 120°. 

cation. The effect is particularly amenable to investigation 
using the inverse ray shooting technique. Figure ^ shows 
that as the source star size increases, the perturbation be- 
comes smaller and broader. Comparison with Figure ^ shows 
that that variations of a planet's mass and projected orbital 
radius can produce qualitatively similar results. This is pos- 
sibly the major shortcoming of the microlensing technique. 
Unless independent information is available to pin down the 
radius of a source star to ~ 20% accuracy, then the precision 
with which the mass and the projected orbital radius of a 
planet can be determined, as given in Section ^, will be sig- 
nificantly compromised. The achromaticity of gravitational 
lensing may be used to advantage in this regard, to obtain 
high quality photometry in several passbands of a source 
star when it is highly magnified. These measurements, com- 
bined with detailed knowledge of reddening along the line 
of sight, should suffice to determine the radius of the source 
star with adequate precision. 

7.2 Source star spots 

The use of microlensing for the detection of small scale struc- 
ture on the source star was recently suggested by Heyrovsky 
& Sasselov (2000). A non-uniform surface brightness distri- 
bution on the source star will produce amplification devia- 
tions from a source star with a uniform brightness distribu- 
tion. These deviations can reach a significant level (> 2%) 
for source-transit events (Heyrovsky & Sasselov, 2000). It 
was noted that fight curve deviations due to source star 
spots may be mistaken for the perturbations due to low 
mass planetary companions of lens star. However Han et al. 
(2000) notes that even if the source star transits the lens, 
the source star spot must pass very close to the lens be- 
fore a detectable deviation will occur. Source star spots will 
therefore affect only a very small number of microlensing 
events. 

In order to confirm that the perturbations due to source 
star spots are indeed negligible compared to those due to 
lens system planets, we performed a number of numerical 



Table 3 . Table of source star spot parameters and maximum light 
curve deviation, max(|<5sp|). These results are for a spot located 
at the centre of the source star disk. 





^max 




max(|<5sp|) 


Rq 


100 


O.lfls 


6 X lO-'^ 


Rq 


100 


O.bRs 


1.3 X 10"* 


2Rq 


200 


O.lRa 


9.1 X 10-5 


2R@ 


200 


0.5iJa 


2 X 10^3 



simulations for typical microlensing geometries. Most mi- 
crolensing events are not source transit events: the minimum 
impact parameter is many times the source star radius. The 
simulations were performed using a range of Mmin, Rs and 
i^spot values, where i?spot is the radius of a source star spot. 
The spot was moved to several locations on the source star 
and the brightness contrast was set at -1-10% with respect to 
the rest of the photosphere. Each light curve using a spotted 
source star was compared to a spotless source star. The devi- 
ation fight curve, 6sp = (Asp__Asl^ computed, where Asp 
and As are the amplification light curves of a spotted and 
spotless source star respectively. For each simulation, the 
maximum value of |5| was recorded. The results are shown 
in Table ^. As we can see from the table, only in the most 
extreme case might the perturbation due to a source star 
spot approach the magnitude due to a low mass planet. 

7.3 Distance and other effects 

In the simulations presented here the parameters Dos, Doh 
and Ml have been set at default values of 8 kpc, 6 kpc and 
0.3Mq respectively unless otherwise stated. These are typ- 
ical values for real events, but they will be subject to large 
event-to-event variations in real cases. A complete analysis 
of the data for actual events will require direct knowledge of 
these parameters. Also, observations of an event during the 
critical time tpwHM only, although sensitive to the presence 
of planets, will not suffice to determine accurately the value 
of j4niax for an event. Knowledge of this quantity is required 
to quantify planetary parameters. 

Additional measurements will be needed in actual 
events to determine the values of the above parameters. 
The distance to the source star Dos may be determined 
approximately from the measurements described above for 
determining its radius. The remaining quantities may be 
evaluated by directly observing the lens and source stars 
a few years after a microlensing event, when they have 
diverged sufficiently to be resolved (Alcock et al, 2001). 
Multi-passband photometry of the lens star would enable 
its type and mass Ml to be determined, and also its dis- 
tance DoL approximately. At the same time, the baseline 
intensity of the source star could be accurately measured, 
and j4max evaluated accurately. These measurements would 
provide absolute values of all microlensing parameters of the 
event. However, the value of the quantity (1 — d) that appears 
in Equation |l| would be quite uncertain, as it involves the 
difference of two similar but only approximately determined 
quantities, Dqs and Dql • This quantity could be better con- 
strained through direct measurement of vt and use of the 
equation vt ~ Re / iE ■ The net result of these measurements 
would be the absolute determination of the mass and pro- 
jected orbital radius of any detected planet, and also the 
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identification of its liost star type. We note, liowever, tfiat 
a telescope witfi capabilities approaching that of the VLT 
interferometer may be required. 

Additional information on the parameters Dos, Dql, 
and Ml may be obtained through parallax measurements 
(Refsdal, 1966; Gould, 1995; Holz & Wald, 1996). These 
would be most useful if carried out from a satellite in a 
distant (solar) orbit. As depicted in Figure ^ the source, 
lens and Earth may then be treated as collinear at the time 
of peak magnification for a high magnification event. The 
amplification as detected by the satellite at this time could 
be determined from the light curve recorded by the satellite. 
This would yield the impact parameter u of the lens at this 
time as seen from the satellite relative to the Einstein ring 
radius of the lens -Re. Also, the baseline b of the satellite 
could be determined from the satellite's orbit. These quan- 
tities would then satisfy the parallax equation: 



u 
b 



Dls 
Dos' 



(10) 



Equation ^ could be used to constrain the parameters 
Dos, DoL and A/l. 



8 ORBITAL MOTION DURING AN EVENT 

Dominik (1998) first studied the effect of orbital motion of 
a binary lens during a microlensing event. It is generally 
held that for most planetary microlensing events, the ef- 
fect of the orbital motion of the planet around the lens star 
will be negligible. This assumption appears particularly safe 
for high magnification events, where the critical observa- 
tion time tpwHM is so short. To confirm this we determined 
typical parameters for which planetary motion during the 
peak of an event would be detectable. A face-on planetary 
system was assumed with a O.SMq lens star and a Jupiter 
mass planet, e = 3.2 x 10~^ orbiting at 0.05-Re. These pa- 
rameters correspond to a "hot Jupiter" orbiting a typical 
lens star. The light curve shows a clear deviation from a 
model with a static planet, see Figure 0. The light curve 
perturbation due to a moving planet is seen to be broader 
(narrower) than that for a static planet when the direction 



Figure 11. Planetary perturbation broadening due to a "hot 
Jupiter" moving in its orbit during an event. The intersection 
of the planet-star axis with the source star track is moving in 
the same sense as the source star, from left to right. The inset 
shows the dynamic microlensing geometry (not to scale). Using 
the simulation parameters of the text, the angle through which 
the planet moves during tpwHM is A9 ~ 10.7°. The initial value 
of e was 45° . 



of the source star is in the same (opposite) as the moving in- 
tersection of the planet-star axis with the source star track. 
The simulation was repeated for a Jupiter mass planet at 
O.SRe, using vt = 20kms~^ and ylmax ~ 10. Similar results 
were obtained, with the broadening effect apparent during 
the whole peak perturbation. 



9 COMPARISON OF THE HIGH 

MAGNIFICATION AND CAUSTIC 
CROSSING TECHNIQUES 

The first planned searches for planets using the microlens- 
ing technique were not based on the high magnification 
technique that has been described here (Albrow et al., 
1996; Pratt et al., 1996). Rather, they were based on the 
well-known "caustic crossing" technique (Mao & Paczynski, 
1991; Gould & Loeb, 1992; Bolatto & Falco, 1994; Bennett & 
Rhie, 1996; Wambsganss, 1997; Bozza, 1999). Caustic cross- 
ing events, in which the light curve undergoes large, rapid 
changes in amplification, can yield a great deal of informa- 
tion about the lens system and the source star. However, 
even if a planet is known to be in the lens system, the prob- 
ability for a caustic crossing to occur is low. Moreover, the 
time of occurrence of the caustic crossing is unpredictable. 
For these reasons, the observational study of planets by the 
caustic crossing technique has proved to be difficult. As a 
practical technique, it is preferable to be able to concentrate 
telescope resources during a short well-defined time period, 
rather than over a longer, indeterminate period. This is es- 
pecially true for ground-based telescopes that are being used 
for several purposes. 

Ultimately, a wide-angle space-borne telescope dedi- 
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Figure 12. Planet positions for wliicli a caustic crossing oc- 
curs within t G [— tEi*E]- The planet mass ratio for these planet 
positions is e = 10~* and the minimum impact parameter is 
Wmin = 0.01. The position axes are in units of R^. Planet posi- 
tions where the source star track crosses the central caustic are 
plotted in grey; those where the track crosses a planetary caustic, 
with dark shading. Lightly shaded regions denote planet positions 
where both a planetary and central caustic crossing occurs. The 
solid line indicates the regions where a planet of the same mass 
can be detected via the high magnification technique, assuming 
a conservative detection accuracy of 1.0%. 

cated to microlensing could achieve the greatest sensitivity, 
by continuously monitoring a very large number of main 
sequence stars. This would yield sensitivity to terrestrial 
planets over a significantly wider range of orbital radii than 
the high magnification technique can achieve. The Galactic 
Exoplanet Survey Telescope (GEST) has been proposed to 
exploit this capability (Bennett & Rhie, 2002). 

In order to further compare the caustic crossing and 
high magnification techniques, we estimated relative fre- 
quencies of planetary detections that may occur with and 
without caustic crossings in events of high magnification. 
From the perturbative analysis of Bozza (1999), the caustic 
curve maps for a range of single planet lens systems were 
generated, varying planet position, mass and the minimum 
impact parameter, Mmin. For each set of planet parameters, 
it was determined if a caustic crossing occurred at any time 
within the interval t £ [— tE,tE]. A solar radius source star 
was used. In order to preserve the perturbative nature of 
the analysis of Bozza, parameter sets where the first order 
perturbative term exceeded 0.05 were ignored. An example 
image of the results obtained in this analysis is shown in 
Figure |l^. The positions where the source star track crosses 
a planetary and/or the central caustic are shown. From this 
it is evident that nearly all planetary detections in high mag- 
nification events will not involve caustic crossings. 



10 CONCLUSION 

The investigations presented here provide further support 
for the practicality of studying planetary systems in gravi- 
tational microlensing events of high magnification, including 
planets with masses as low as that of Earth. 

It has been shown that the perturbations produced by 
planets with position angles 30° < < 150° or 210° < 
6 < 330° with respect to the source star track occur in 
the restricted time interval [— 0.5fFWHM, O.StFWHM] about 
the peak of the hght curve. Here, fpwHM = -j^^, and is 
typically ~ 14 hours when Amax ~ 100. This implies that the 
planetary deviations of two thirds of all planets occur in this 
relatively restricted time interval, a significant advantage for 
carrying out practical observations. 

It has also been demonstrated that the position angle 
of a planet is simply related to the time and sign of the 
perturbation it produces on the fight curve, and that its 
mass and projected orbit radius may be determined from the 
height and shape of the planetary deviation, provided the 
radius of the source star can be determined independently. 

It has been shown that Earth-mass planets are typi- 
cally detectable with 1-m class telescopes provided they lie 
within about 25% of the Einstein ring radius of the lens 
star, i.e. from about 1.5 - 2.5 AU. Neptune-mass are typ- 
ically detectable provided they lie beyond about 0.8 AU, 
and Jupiter-mass planets are detectable almost anywhere. 
These zones of detectability, whilst large, only marginally 
include planets in the "habitable zone". It was also shown 
that nearly all planetary detections in high magnification 
events will not involve caustic crossings. 

Multi-planet systems can be systematically charac- 
terised. To a good approximation, the joint perturbation is 
the sum of the individual perturbations of the component 
planets. This approximation may be used as a first approx- 
imation in planetary modelling. 

Source star effects have been investigated. The main ef- 
fect is the size of the source star. This may be able to be 
determined by multi-passband photometry carried out dur- 
ing the course of a microlensing event. The effects of spots 
on the source star are likely to be negligible. The effect of 
orbital motion of a planet during a microlensing event has 
also been investigated. In most cases this is entirely negli- 
gible, although in the case of hot Jupiters it is likely to be 
significant. 

Follow-up observations carried out a few years after a 
microlensing event when the lens and source stars are re- 
solvable should suffice to characterize the parameters of any 
detected planet, and its host star, absolutely. 

Finally, we have found that the numerical inverse ray 
shooting technique for analysing microlensing events can be 
conveniently carried out on a cluster computer, and that 
complicated physical situations can thereby be rapidly sim- 
ulated. 
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APPENDIX A: INVERSE RAY SHOOTING 
CODE DETAILS 

An improved version of the inverse ray shooting code has 
been developed which utihses a special property of the mi- 
crolensing mapping functions. A rectangular region Bm is 
mapped to a new region B!^ in the source plane, see Figure 



Al 



It can be shown that for any lensplane point L — {x, y) 
contained in a lensplane boundary box Bm , the correspond- 
ing source plane point S = {^{x,y),ri{x,y)) will be con- 
tained within the boundary B'^ which is the smallest rect- 
angle in the source plane containing B'^ ■ The proof of this 
statement follows by demonstrating that, for regions Bm 
not containing any lenses, the mapping functions ^ and ri do 
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Figure Al. Boundary box theorem: Any lens plane point L con- 
tained within boundarjy^ box Bm , will be mapped to a point S 
contained within box B' . 



not have any local extrema. If none of the elements of B'm 
are contained within the source star track region of interest, 
the area of the lensplane defined by Bm can be discarded. 
The improved version of the code has a run rate almost one 
order of magnitude less than the original version. This is 
significant when modelling actual datasets. 



Microlensing mapping functions 

Due to the gravitational microlensing, any given point 
L(x,y) from the lens plane is mapped into a single point 
T)) in the source plane: 



nix, y)^y- J2k^i ^ 

k 

where: 

Axfe ^ X — Xk = Tk cos ipk, 
^Vk = y - yk = r-k sin ipk , 



(Al) 



(A2) 



and Tk and ipk are the corresponding polar coordinates. The 
mapping functions f (a;, y) and rf^x, y) are defined everywhere 
in the lens plane, except at the lensing points themselves. 
We shall exclude these singular points and consider only the 

case Vk 7^ 0. 

The microlensing mapping ( |Al[ ) can also be regarded as 
a single complex function: 



JV 



2fe 



= 1 



Zk 



N 

fc=i 



nik 



Zk 



(A3) 



Clearly, (^{z) is not an analytic function (i.e. it is not 
complex-diff'erentiable) , since it depends on both z and z. 
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However, one can always differentiate this function partially 
with respect to x and y (or even z and z). 

Partial derivatives 



Starting from Equation (A3) it is easy to calculate the first 
partial derivatives of ({z) with respect to z and z: 



dz ' 



az ^-^ 



(A4) 



Using the relative polar coordinates Vk and (pk, the first 
derivatives can be expressed as: 



dz 
where 



= C2+iS2, 



N 

\ - cos nifk 
Cn = y rrik — , 



N 

\ - sm rupk 

fc=l 



(A5) 



(A6) 



From Equations (A5) one can also calculate the derivatives 
with respect to x and y: 

| = (1 + C.)+.S„ 



1^ = S2 + i(l-C2). (A7) 
oy 



Higher derivatives of ^(2) can be obtained by further 
differentiating the first derivatives (A4). Since d(/dz — 1, 
all higher derivatives with respect to z are zero. Any mixed 
derivatives with respect to both z and z are equal to zero as 
well, since d(/dz does not depend on z. The only remaining 
case involves higher derivatives with respect to z alone: 



(A8) 



^^(_l)»+isiV — 

dz' ^ ' {z-ZkY+'-' 

fc=i 

This can be rewritten in terms of the relative polar coordi- 
nates Tfc and (fk as: 

9X 



dz 



(A9) 



One can also obtain the following expression for higher 
derivatives with respect to x and y: 



(AlO) 



dxPdyi dz'' 
where s — p + q and s > 1. 

Stationary points of the mapping functions 

In order to determine whether the microlensing mapping 
functions ^(a;, y) and ri{x, y) have any local extrema, one has 
to locate the stationary points, where the first derivatives are 
equal to zer o. B y extracting the real and imaginary parts of 
Equations (A7), one obtains: 



9L 
9y 



S2 



dr] 
dy 



C2. 



(All) 



The stationary points are found by simultaneously solving 
S2 ~ Q and C2 = — 1 for ^{x,y), and 52 = and C2 = 1 
for r){x,y). Obviously, the two functions cannot have any 
common stationary points. 



It is well known that the nature of a stationary point of 
a real function f{x,y) can be determined by examining the 
sign of the Hessian determinant: 



D 



dx'^ dy^ 



dxdy 



(A12) 



Using Equation (AlO) for s — 2 and then separating the 
real and imaginary parts, one can demonstrate that both 
functions ^{x,y) and rj{x,y) give the same value for D: 



D^-A{Cl + Sl). 



(A13) 



This excludes any possibility of D being positive, but one 
might still have -D = 0, if all second derivatives are equal 
to zero. This means that further examination is needed to 
determine the type of a stationary point. 



Taylor's formula 

One can use Taylor's formula to express the values of the 
mapping functions at a given point L{x,y) with respect to 
a nearby fixed point Lo{xo, yo)'- 



m = e(io) + e::i 

where the differential operator d" is defined as: 



d =[Ax-+Ay- 



(A14) 



(A15) 



The reminders d"^ and d"r] are evaluated at two different 
points, G and H , located somewhere along a section of the 
straight line connecting Lo and L. 

Let us assume that Lo is a stationary point for one of 
the mapping functions. Since the first partial derivatives are 
zero at Lq, the summation in Taylor's formula can effectively 
start from s = 2. In addition, we are only interested in the 
'problematic' case when all second derivatives are also equal 
to zero. This means that the summation can actually start 
from s = 3. In general, we have to allow a possibility that 
all derivatives up to the order (n — 1) at Lo are equal to 
zero. However, there must exist a high enough value for n so 
that at least one of the n-th partial derivatives is not zero. 
This means that either C„+i{Lo) or Sn+i{Lo) is not zero. 
Taylor's formula can be rewritten in the following form, by 
keeping the reminders only: 

= ?(i) - C(io) = = ^ 3le<5"C(G), 

Arj = viL) - v{Lo) = = i 3m S-CiH), 

where 5"^ is the complex differential operator: 
d 



(A16) 



/ A o ^ d\' 



s" = 

It can be shown that: 

— jn] {Cn+i + iSn+i) (cosn?/) — isinn?/;). 



(A17) 



(A18) 



where 7 — (— l)""'"'^p", p and ip are the polar coordinates of 
Az and C„+i and S„+i are evaluated either at G (for A^), 
or at H (for Arj). By extracting the real and imaginary 
parts of 5"^, one can obtain the following expressions for 
the residuals and Arj: 
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= 7 TCn+i (G) cos + S„+-i (G) sin nt/;] , 1 

(A19) 

Ar; = 7 [S^+i [H) cos - C„+i [H) sin nV"] . J 

If a mapping function has a local extremum at io, then 
the corresponding residual must have the same constant 
sign at every point L around Lo. It is easy to demonstrate 
that neither A^ nor A?) can keep the same sign everywhere 
around Lo. For example, let us consider A^. We know that 
at least one of the coefficients C„+i and Sn+i is not zero 
at Lo- Let C„+i(io) 7^ 0. One can always pick a point L\ 
at i/ii = 0, so that A^i = 7C„+i(Gi). Another point, L2, 
can be at ijj2 = 180°/n, so that A^2 = — 7C„+i(G2). We 
can always choose the radius p to be small enough, so that 
C„+i(Gi) and C„+i(G2) are of the same sign depending 
on the sign of C„+i(Lo), since C„+i is continuous at Lo- 
in such a case, we have opposite signs for A^i and A^2, so 
that Lo cannot be a local extremum. If, on the other hand, 
we start with Sn+i{Lo) 7^ 0, then L\ and L2 can be set 
at Vi = 90° /n and '4>2 = 270° /n, respectively, so that the 
cosine function is zero and Afi and A^2 have opposite signs 
again, for the same sign of Sn+i{Gi) and Sn+i(G2)- 

Using the same idea, it is easy to demonstrate that Ar; 
cannot be of the same sign everywhere around Lo either. 
Therefore, the microlensing mapping functions cannot have 
any local extrema. 
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